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Abstract 

The generic unparticle propagator may be modified in two ways. Breaking the 
conformal symmetry effectively adds a mass term to the propagator, while con- 
sidering vacuum polarization corrections adds a width-like term. Both of these 
modifications result naturally from the coupling of the unparticle to standard 
model (SM) fields. We explore how these modifications to the propagator af- 
fect the calculation of the lepton anomalous magnetic moment using an integral 
approximation of the propagator that is accurate for d < 1.5, where d is the 
unparticle dimension. We find that for this range of d and various values of the 
conformal breaking scale /z, the value of g — 2 calculated when allowing various 
SM fermions to run in the unparticle self-energy loops does not significantly de- 
viate from the value of g — 2 when the width term is ignored. We also investigate 
the limits on a characteristic mass scale for the unparticle sector as a function 
of jj, and d. 



1. Introduction 

Unparticle physics [l|, Q represents an intriguing possibility for physics be- 
yond the standard model. In this scenario, there is a new sector which is 
conformal below a mass scale Ay. It can be described by an effective field 
theory of so-called "unparticle" operators, which in some respects imitate a 
non-integer number of massless particles. The phenomenology of unparticles 
has been studied in great detail. In particular, the anomalous magnetic mo- 
ments of the electron and the muon have been calculated in various unparticle 
scenarios 0, 13, 0, 0] • Here we explore the effect of adding a width-like term 
to the denominator of the propagator on the calculation of electron and muon 
g — 2. We also study in detail the relationship between the conformal breaking 
scale (i and these quantities. 

The propagator for unparticles can be obtained from the spectral represen- 
tation 

//■oo t 7\/r2 
d\e^(O u (x)Ol(0)) = j Q — pu (M 2 ) p2 _ M2+ . e , (1) 
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Figure 1: The imaginary part (left panel) and real part (right panel) of the 
unparticle propagator Au(p 2 ) for various values of unparticle dimension d, in 
the absence of unparticle sector conformal symmetry breaking. 



where Ou is an unparticle operator and pu(M 2 ) is the unparticle spectral func- 
tion. If Ou has scaling dimension d, conformal invariance determines 



P u{M z ) = A d {M A f- z , (2) 

where 



A d-77^V2d— 7vn7o3T ( 3 > 



167T 5 / 2 T(d + 1/2) 

(^pr(d- i)r(2d) 

is chosen so that the rf — > 1 limit gives the propagator for a normal particle. Fol- 
lowing Q we will take d S (1,2). Inserting this spectral function into Eq. Q]and 
evaluating the integral yields the propagator 

PU(P 2 ) = T ^ T A-P 2 - I er-*. (4) 
2 sin(7ra) 

The above is valid for a scalar unparticle. For vector unparticles, which we 
consider in this paper, this propagator must be multiplied by the projector 
— {diiv + bp^p v /p 2 ). The constant b is in principle undetermined because the 
unparticle theory need not have a gauge symmetry, but it will not affect the re- 
sults of any of our calculations so we set it to zero from now on. The propagator 
is then P^(p 2 ) = -g^Puip 2 ). 

To get a feel for the behavior of this unusual propagator, we show the imag- 
inary and real parts of Au(p 2 ) = —iPu(p 2 ) for positive p 2 in Figs. [Tal and [TBI 
respectively. It will be convenient for us to use Au rather than Pjj in our anal- 
ysis and we will often refer to Au as "the propagator." It should be noted that 
the imaginary part of Au{p 2 ) is zero for negative p 2 , while the real part is not. 

Following [8| we shall write the coupling of vector unparticles to Standard 
Model (SM) fermions in the form 

^ = *-^7m«- (5) 

A 3 
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The parameter A3 is not fundamental; it can be expressed as the combination 



where M and ciyy are parameters of the underlying UV theory as described 
in Q. Since the coupling c and the scale A3 only appear in the combination 
given in Eq. and we have no theoretical guidance as to what c may be (other 
than naturalness, which might suggest that c is 0(1)), in our analysis we will 
take c = 1 and let the coupling strength vary by varying A3. 

An interesting modification of the unparticle scenario involves breaking the 
conformal symmetry of the unparticle sector at a scale /1. This can result from 
interactions between the unparticle and the Higgs, for example. An ansatz for 
the spectral function in this case essentially consists in removing modes below 
p, 2 from the spectrum, yielding Q 

H {M 2 ) = A d 9(M 2 - fj, 2 )(M 2 - f i 2 ) d - 2 . (7) 

With this spectral function, the propagator becomes 

For the fractional power, we take the standard branch cut such that 

In Fig. [3al we show the imaginary part of A^(p 2 ) = — iP^(p 2 ) for d = 1.3 and 
different choices of //. It rises steeply near p 2 = fj, 2 and vanishes for p 2 < fi 2 . 
In Fig. [2b] we display ReA^(p 2 ) for fi = 4 GeV for different values of d. Unlike 
the imaginary part, the real part is nonzero for p 2 < /i 2 , though for d = 1.5 it 
is zero for p 2 > /i 2 . 

Some of the strongest experimental constraints on unparticle physics can be 
evaded with a sufficiently large value of \x @, Ho| . To understand this, note 
that when there is a nonzero fi, the form of the spectral function (see Eq. [7J 
is such that low M 2 modes no longer contribute to the propagator (cf . Eq. Q]) . 
It is well-known that the bulk of the contribution to g — 2 in QED comes 
from low M 2 modes in the propagator 

El Ei El; hence . for lar s 

e enough 



values of /z, the contribution to g — 2 from unparticles is negligible. It has been 
suggested JlOj] that the propagator given in Eq. [8] should be further modified 
by summing the contributions to the propagator from diagrams containing an 
arbitrary number of vacuum polarization bubbles, as is done when finding the 
renormalized mass and width of a "normal" particle. It will be interesting to see 
how this modification to the propagator affects the above intuition regarding 
the size of g — 2. In particular, we wish to know whether such a modification, 
which should lead to a nonzero spectral function for values of M 2 less than /x 2 , 
will result in a significantly larger value of g — 2. 
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-ImA(p ) with broken scale invariance, d=l. 
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ReA(/? ) with broken scale invariance, d=1.3 
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Figure 2: The negative imaginary part (left panel) and real part (right panel) of 
the unparticle propagator with broken conformal invariance A^(p 2 ) for d = 1.3. 
The negative imaginary part is shown for four different values of the conformal 
breaking scale /i, while the real part is shown for five different values of the 
unparticle dimension d. 



In this paper we calculate the contribution to muon g— 2 from a vector unpar- 
ticle both with and without vacuum polarization corrections to the propagator. 
We obtain bounds on d, /i, and A 3 from the experimental limits on muon and 
electron g — 2, and investigate how these limits are affected by these corrections. 
Specifically we will consider the modification to the unparticle propagator from 
vacuum polarization graphs with fermions running in the loop. In particular, 
we will consider the following 3 cases: 

1. The unparticle couples only to the lepton for which we are calculating 
.9-2. 

2. The unparticle couples to all leptons with equal strength, but to no other 
particles. 

3. The unparticle couples to all standard model fermions with equal strength, 
but to no other particles. 

The calculation of the unparticle contribution to g — 2 in these cases will re- 
quire obtaining an approximate dispersion integral representation of the modi- 
fied propagator, whose validity and accuracy we will investigate. 



2. Vacuum polarization corrections to the propagator 

We will first consider the modification of the propagator from vacuum po- 
larization diagrams. While in principle we could consider decays to scalar or 
vector particles, and possibly even other species of unparticle, for simplicity we 
confine our attention to loops involving standard model fermions. For a given 
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fcrmion in the loop, the amplitude for the vacuum polarization graph is 




(10) 



where 



(ii) 



and B(f3) is the Feynman parameter integral 



B{f3)= [ dxx(l-x)log(l/(l-4Px(l-x))) 



(12) 



Jo 



Our unparticle theory is not a gauge theory. Since the rest of the g-2 calculation 
involves only QED, however, the contribution proportional to q^q v will vanish 
by the Ward identity. 

Summing over diagrams with an arbitrary number of vacuum polarization 
loop insertions yields the propagator P!t u {p 2 ) — — P;ji(p 2 ) , with 



To see the effect of this modification to the propagator, in Fig. [3a] we show 
-ImA^(p 2 ) = -luii-iPjjfip 2 )) for d = 1.1, n = 10 GeV, and A 3 = 1 TeV 
with four different choices of vacuum polarization corrections. The peak near 
p 2 = /i 2 is broadened considerably by the corrections, showing that the AdJl(p 2 ) 
term in the denominator of the propagator has an effect similar to that of a 
decay width for a normal particle. In Fig. I3b[ on the other hand, we show 
— ImA^(p 2 ) again for d = 1.1 and fx = 10 GeV, but now we only consider 
vacuum polarization corrections from all SM fermions and let A3 take on four 
different values. Because decreasing A 3 effectively increases the unparticle-SM 
coupling, the propagator peak is broadened more for lower A 3 . 

3. Calculating g — 2 using a dispersion integral 

Now that we have the form for the unparticle propagator with vacuum po- 
larization corrections, we turn our attention to the question of how to calculate 
the one-loop contribution to g — 2 from unparticles both with and without these 
corrections to the propagator. This contribution comes from the diagram in 
Fig.Hl The amplitude iM? corresponding to this diagram satisfies 



iA d 



(13) 



2sm(7rd)(^ 2 - p 2 - ie) 2 - d - A d U(p 2 ) ' 



/ 



d A k uip')^^' + mV - 2m(k + k'Y)u(p) 
(2tt) 4 (fc 2 -m 2 + ie){k' 2 - m? + ie)D u 



(14) 
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-ImA(/5 ) for d— 1.1, A3-I TeV -ImA(/? 2 ) for d= 1.1 with vacuum corrections from all SM fermioi 
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Figure 3: The negative imaginary part of the propagator, for d — 1.1 and \x = 
10 GeV, for A3 = 1 TeV and different choices for loop corrections (left panel) 
and for various values of A3 (right panel) . In the left panel, the "uncorrected" 
curve shows the propagator in the absence of vacuum polarization corrections, 
while the other three curves show the value of the propagator if only muons, 
all leptons, or all SM fermions are allowed to run in the loops of the self-energy 
terms which lead to a width-like term in the unparticle propagator. In the right 
panel, different curves represent different values of A3; since this parametrizes 
the fermion-unparticle coupling, it affects the shape as well as the normalization 
of the propagator if the propagator includes the width-like factor from vacuum 
polarization corrections. 




Figure 4: The unparticle correction to the muon g — 2. The photon momentum 
q is incoming and the muon momenta follow the fermion lines. 
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where Du is given by 



D u = { l ?-(p-kf-ief- d (15) 
for an unparticle propagator without vacuum polarization corrections, and 

D U = (p 2 (p kf le f- d —^—U(( P kf) (16) 

for an unparticle propagator with vacuum polarization corrections. 

Note that the numerator is identical to that which one obtains in the calcu- 
lation of the QED one-loop contribution to g — 2. Since we are only considering 
vector unparticles, which conserve parity, we can write the amplitude as 

iM» = ieu{p') (V*i(<7 2 ) + i^F 2 (q 2 )) u(p) . (17) 
\ 2m J 

The unparticle contribution to the quantity a = q/2 — 1 is then given by a u — 
F 2 (0). 

When calculating the QED one-loop contribution to g — 2 (where the particle 
in the loop is a photon rather than an unparticle, and hence Djj — > (p — k) 2 + ie), 
we can use the method of Feynman parameters to combine the three factors in 
D and shift the integration variable to I = k + yq — zp. The result is 

9 9 f d 4 l f 1 z(l-z) 
F 2 (0) = 8*Wy ^ I d X dy dZ (l2 } A + > e)3 . (18) 

Note that there are no powers of I in the numerator of the integrand. As a 
result, the integral is convergent. The simplifications (such as determining that 
terms proportional to I contribute to Ft rather than F2) that make this explicit 
require the use of Feynman parameters. 

In the case of unparticles without vacuum polarization corrections to the 
propagator, the Feynman parameter method can still be used, though in this 
case the denominator is raised to the power 4 — d. With the corrected propaga- 
tor, however, it is impossible to combine the terms in the denominator into an 
expression quadratic in momentum variables. Hence, we cannot use Feynman 
parameters in this case to obtain a relatively simple, manifestly convergent inte- 
gral expression for F 2 (0) analogous to Eq. [TSl It is difficult to imagine that one 
can calculate the integral in Eq. [14] analytically in this case. One could hope to 
perform the integral numerically, but without the simplifications obtained using 
Feynman parameters, the integral has divergent contributions whose cancella- 
tion would be difficult to guarantee in a numerical calculation. We therefore 
seek to recast the integral in a form amenable to the use of Feynman parameters. 

To do this, we try writing the unparticle propagator as a dispersion integral 
(cf. Eq. [T]) which contains the standard vector boson propagator. Because 
p(M 2 ) is independent of the loop momentum, it can be pulled out of the loop 
momentum integral. The vector boson propagator, however, remains under the 
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loop integral and we end up with 



f°° rIM 2 

a U = I ^p(M 2 )I(M% (19) 

where p(M 2 ) is the spectral function and I(M 2 ) is the contribution to g — 2 
at one loop from a vector boson with mass M and coupling to fermions given 
by Eq. [5] Specifically we have that I(M 2 ) is given by the Feynman parameter 
integral 

"1 x 2 (1 ^2 m 2 



I{M ) = i^Ap* Jo d W+(l-x)M3 ' (20) 

where m is the mass of the fermion for which we are calculating g — 2 (generally 
the muon), and as noted above we will set c=l in what follows. The Feynman 
parameter integral in Eq.f^Olcan be performed to obtain an analytical expression 
for I(M 2 ), leaving the M 2 integral in Eq. [12] to be done numerically. 

This method of expressing the propagator as a dispersion integral clearly 
works in the case without vacuum polarization corrections, since we already 
have expressions given above for p{M 2 ) (Eqs.[2]and[7|). Of course in this case, 
as mentioned above, one could do the calculation using Feynman parameters 
0, [|| without recourse to the dispersion integral. This provides a useful check 
on our results. 

To use this method when we include vacuum polarization corrections to the 
unparticle propagator, however, we still need to express that propagator as a 
dispersion integral. In the next section, we will study whether this can be done. 



4. Expressing the propagator as a dispersion integral 

We need to find an expression for the spectral function of the vacuum 
polarization-corrected unparticle propagator. To this end, we first note that the 
spectral function for the "tree-level" propagator in both the unbroken (Eq. \2§ 
and broken (Eq. [7]) conformal symmetry cases satisfies the relation p(M 2 ) = 
2ImA(Af 2 ). In other words, the spectral function is (twice) the imaginary part 
of the propagator. This is actually due to a general result in complex analysis, 
which states that for any function of a complex variable f(z), 

m . -i r jz&u, , pi, 

n J-oo z-z' + «e 

if / is analytic in the upper half-plane. (The difference between the — factor 
in this equation and the -i- in Eq. [T] accounts for the factor of 2 in the relation 
p(M 2 ) = 2ImA(M 2 ), while the i in the numerator of Eq. [1] combined with the 
i in A = ImP lead to the lack of a relative sign.) The uncorrected unparticle 
propagator is indeed analytic in the upper half-plane in both the broken and 
unbroken conformal invariance cases. 

It therefore seems like a simple matter to obtain the spectral function for 
our vacuum polarization corrected unparticle propagator, Eq. 1131 It should be 
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Pjjj{M 2 ) = 2ImA^(M 2 ). The problem is that it is not clear if the corrected 
propagator A^(M 2 ) is analytic in the upper half-plane. If it is not, this relation 
may not hold. 

We would like to test if this condition holds exactly, indicating that the 
propagator is analytic in the upper half-plane. It may also be possible that 
the propagator is not analytic in the upper half-plane, but is close enough to 
an analytic function that the dispersion integral with p^(M 2 ) = 2ImA^(M 2 ) 
serves as a good approximation to the propagator. To investigate these possi- 
bilities, we adopt a numerical strategy. We will assume pu(M 2 ) = 2ImA^(M 2 ) 
as a putative definition for the spectral function associated with the vacuum 
polarization-corrected unparticle propagator. We then want to know if the dis- 
persion relation 

/•OC JTijfl 1 

2W)I-jf < 22 > 

holds. Note that the lower range of integration may be taken to be because 
ImA^(p 2 ) is zero for p 2 < 0. The left-hand side of this equation is the true value 
of the propagator, which we can simply evaluate using Eq. [131 The right hand 
side is our putative spectral representation for the propagator, and we can com- 
pute its value by carrying out the integral numerically. We can then compare the 
two sides. If they do not agree, we learn that, for the loop-corrected propagator, 
the relation Ptf(p 2 ) = 2ImA^(p 2 ) does not in fact hold (or, equivalently, that 
Eq. [21] does not hold) , because the propagator is not analytic in the upper- half 
plane. If the two sides of Eq. [22] do agree to a high precision, we can conclude 
that the propagator is analytic (or at least well-enough approximated by an 
analytic function) in the upper half-plane, so that the RHS provides the correct 
spectral representation (or at least good approximation) of the true propagator. 
If this is the case, we can use p^{p 2 ) as a good approximation to the spectral 
function in the g — 2 calculation. 

It should be noted that, following [HI, when the mass m of the fermion in 
the loop satisfies 2m > p, there will be a value of M 2 for which the denominator 
of A^(Af 2 ) is zero. This leads to an additional delta function term in p^{M 2 ). 
Such a value can only occur when the loop contribution II(M 2 ) and the first 
term in the denominator of the propagator, (2sm{ird) / Ad)(p 2 — p 2 — ie) 2 ~ d , are 
both real. The loop contribution from a given fermion with mass m is real for 
M 2 < 4m 2 while the first term in the denominator of the propagator is real for 
M 2 < p 2 . Hence, the denominator cannot vanish for p 2 > M 2 > 4m 2 , where 
m is the mass of the lightest fermion running in the loop, unless there were 
some fortuitous cancellation of the imaginary parts of the loop contributions 
from the various fermion species. It can be easily verified that there is no such 
cancellation. It may be further shown that, given the signs of the imaginary 
parts of the various terms in the denominator of the propagator, it also cannot 
have a zero for M 2 > p 2 . For M 2 < Am 2 , where m is the mass of the lightest 
fermion running in the loop, both terms in the denominator are strictly real. 
In the situations we consider, however, where the lightest fermion in the loop is 
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either an electron or a muon, and p, takes values of 10 or 100 GeV, the absolute 
value of the first term in the propagator, in this region, is always much larger 
than the absolute value of the loop correction term, so no cancellation can occur. 
Thus, in the situations we consider, there are no delta function contributions to 
the unparticle spectral function. 

In Table [H we compare the actual loop-corrected propagator (Eq. fT5|) to its 
putative spectral representation (RHS of Eq. [22f for a large number of choices 
for d, /it, and p 2 . For each choice of these parameters, we show the fractional 
error between the spectral integral approximation for the pro pag ator and its 
true value. That error is obtained by first, using Mathematica [161 ], performing 
the numerical integration to evaluate the RHS of Eq. [5U (which we'll denote 
by A° lsp (p 2 )). We then compare this to the actual value of the propagator, 
obtaining the fractional error 



rdisp 



A dis P^2 



{p 2 )-^{p 2 ) 



(23) 



This fractional error is obtained for four cases. First, for the case of no loop 
corrections to the propagator. In this case, we know the spectral representation 
with Pj/i{p 2 ) = 2ImA^(p 2 ) to be exact, so we expect zero error, and indeed find 
errors small enough to be attributed to inaccuracy in the numerical integration 
(except when the value of the propagator is exceedingly small, as explained 
in the caption). In the last three columns, the fractional error is shown with 
three different sets of fermions allowed to run in the loop in the unparticle 
propagator — muon only, leptons only, and all SM fermions. We notice that for 
all these cases and for all values of d < 1.5, the fractional error is at most a few 
percent. Thus for these values of d, we can use Ptf(p 2 ) = 2ImA^(p 2 ) to carry 
out the calculation of the unparticle contribution to g — 2 with confidence that 
the result will be accurate. 



5. Results 

We now have a dispersion representation for the unparticle propagator. It is 
exact in the cases with unbroken and broken conformal symmetry but without 
vacuum polarization corrections. In the case with vacuum polarization correc- 
tions, it is a good approximation for d < 1.5. We can then numerically integrate 
Eq. [19] to obtain the unparticle contribution to g — 2. 

Doing this, we find, in fact, that the effect of the propagator corrections on 
the g — 2 contribution is small. Table[2]shows the relative difference between the 
unparticle contribution to g — 2 with and without the inclusion of vacuum polar- 
ization corrections to the unparticle propagator. The difference is negligible for 
all parameter choices for which the dispersion representation of the propagator 
is accurate. Therefore, in the rest of this section, we will discuss results using 
the uncorrected unparticle propagator, knowing they will still hold when the 
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d, /i, p z 


Value 


None 


Muon 


Leptons 


All 


1.2,10,5 


-1.2 x 10"^ 


-8. x 10~ 14 


-1.3 x 10~ 8 


-1.2 x 10~ Y 


-1.7 x 10~ b 


1.2, 10,500 


3.1 x 10~ d 


6.6 x 10" io 


4. x 10~ a 


5.4 x 10~ Y 


6.7 x 10~ b 


1.2, 10,50000 


6.5 x 10~ b 


-1.1 x 10" la 


1.9 x 10" ' 


2.2 x 10~ b 


3.2 x 10~ J 


1.2, 100,5 


-2.9 x 10~ 4 


6.5 x 10" 12 


-3.3 x 10~ 8 


-2.1 x 10~ b 


-2.1 x 10~ 4 


1.2,100,500 


-3. x 10~ 4 


6.6 x 10~ X2 


-2.1 x 10~ s 


-2. x 10~ b 


-2.1 x 10~ 4 


1.2, 100, 50000 


7.8 x 10~ b 


-1.1 x 10" la 


8.3 x 10~ K 


7.9 x 10~ b 


8.4 x 10~ 4 


1.4, 10,5 


-1.4 x 10^ 


-1.3 x 10~ xb 


-7.3 x 10~ b 


-4.9 x 10~ b 


-3.8 x 10~ 4 


1.4, 10,500 


1.9 x 10 - ' 1 


4.6 x 10" ia 


5.6 x 10~ b 


3.7 x 10~ 4 


2.9 x 10" a 


1.4, 10,50000 


1. x 10~ 4 


-3.4 x 10" ia 


1. x 10~ a 


6.8 x 10~ a 


5.6 x 10~ 2 


1.4, 100,5 


-8.8 x 10~ 4 


4.2 x 10" ia 


-8.2 x 10~ b 


-4.9 x 10~ 4 


-3.4 x 10~ 3 


1.4, 100, 500 


-9.1 x 10~ 4 


4.3 x 10~ xd 


-8. x 10~ b 


-4.8 x 10~ 4 


-3.3 x 10" a 


1.4, 100, 50000 


1.2 x 10~ 4 


-2.1 x 10" ia 


6.1 x 10~ 4 


3.7 x 10~ a 


2.6 x 10~ 2 


1.6, 10,5 


-1.9 x 10^ 


9. x 10~ xb 


-9.7 x 10~ 4 


-2.3 x 10~ a 


-6. x 10~ a 


1.6, 10,500 


-3.3 x 10~ ;l 


-4.6 x 10" ia 


-5.6 x 10~ ;l 


-1.3 x 10~ 2 


-3.4 x 10~ 2 


1.6, 10,50000 


-4.9 x 10~ 4 


3.4 x 10" 13 


-3.8 x 10~ 2 


-9. x 10~ 2 


-2.3 x 10~ x 


1.6, 100,5 


-3. x 10" a 


-4.4 x 10" ' 


-5. x 10~ J 


-1.1 x 10~ 2 


-2.8 x 10~ 2 


1.6, 100,500 


-3.1 x 10~ d 


-4.6 x 10"' 


-4.9 x 10 - ' 1 


-1.1 x 10~ 2 


-2.7 x 10~ 2 


1.6, 100, 50000 


-5.3 x 10~ 4 


2.1 x 10" la 


-2.8 x lO"^ 


-6.3 x 10~ 2 


-1.6 x 10~ x 


1.8, 10,5 


-3.5 x 10"^ 


-2.7 x 10" ' 


-4.1 x lO"^ 


-5.7 x 10"^ 


-8.2 x 10~ 2 


1.8, 10,500 


-2.1 x 10"^ 


-1.1 x 10~ x;i 


-6.8 x lO"^ 


-9.4 x 10~ 2 


-1.3 x 10~ x 


1.8, 10,50000 


-8. x 10" a 


8.1 x 10~ X4 


-1.8 x 10~ x 


-2.5 x 10~ x 


-3.5 x 10~ x 


1.8, 100,5 


-1.4 x 10"^ 


-1.1 x 10~ 8 


-9.4 x lO"^ 


-1.2 x 10~ x 


-1.8 x 10~ x 


1.8, 100,500 


-1.4 x lO"^ 


-1.2 x 10~ s 


-9.3 x 10^ 


-1.2 x 10~ x 


-1.7 x 10~ x 


1.8, 100, 50000 


-8.3 x 10~ d 


4.9 x 10~ X4 


-1.5 x 10~ x 


-2. x 10~ x 


-2.9 x 10~ x 



Table 1: The fractional error (5 dls P (see text) in the dispersion integral ap- 
proximation for the unparticle propagator. An entry in the first column gives 
the choice of parameters used for that row, with fj, given in GeV and p 2 in 
GeV 2 . The second column gives the actual value of the propagator (without 
vacuum polarization corrections) for that choice of parameters. The remaining 
columns give the fractional error for the dispersion integral approximation for 
the propagator A^ lsp (p 2 ) using different choices for vacuum polarization correc- 
tions. "None" means no vacuum polarization corrections, "Muon" means only 
the muon loop is included, "Leptons" means the vacuum polarization correc- 
tions from all charged SM leptons are included, and "All" means the corrections 
from all charged SM fcrmions are included. 
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vacuum polarization corrections to the propagator are included^ It is possi- 
ble, of course, that H(p 2 ) is significantly modified at high p 2 by heavy particles 
whose presence we have not taken into account, for example from other hidden 
sectors to which the unparticle is coupled. We note that this possibility could 
cause a larger effect on g — 2 and other observables from vacuum polarization 
corrections. 

Strong experimental constraints exist on g — 2 for both muons and electrons. 
In the case of muons, the measured value of a is three standard deviations larger 



than the calculated SM value. Specifically [17 1 



a^ xp - af 1 = 2.95 ± 0.81 x 1(T 9 . (24) 

As a very conservative upper bound, we require that the unparticle contribution, 
which is always positive in the case of vector unparticles, be no more than three 
standard deviations above the central value, i.e. 

af t < 5.38 x 1CT 9 . (25) 

In the case of electrons, the measured and calculated values are in good 
agreement. We will take the conservative bound [6| a% < 1.5 x 10~ n , which, we 
note, is over five times the error on the most recent experimental measurement 
0. 

In Fig. I5a( we display the unparticle contribution to muon g — 2 in the 
case of broken conformal symmetry as a function of \i and compare it to the 
experimental bounds. As described above, this contribution is obtained by 
numerically integrating Eq. 1191 again using Mathematica. The contribution 
begins to drop quickly as \x is increased above m^, and the statement that for 
\i > 1 GeV, the experimental bound is evaded, is shown to be true for, e.g., 
A 3 = 1 TeV and d = 1.5. It is interesting to note the large hierarchy between 
the d = 1.1, 1.5, and 1.9 contributions. 

Finally, in Fig. I6al we show the lower bounds on d and fj, from muon g — 2 
for various choices of A3. We note the very strong dependence on d and p. In 
Fig. I6bl we display the same bounds but from electron g — 2. It is qualitatively 
similar though the bounds are much weaker. 



6. Conclusions 

In this paper we have examined the correction to lepton g — 2 in several 
unparticle scenarios. In particular, we studied the model in which in addition 
to breaking the conformal symmetry of the unparticle sector at a scale /i, we 
also added a width- like term from vacuum polarization to the propagator. To 



1 In fact we are only certain this is true for d < 1.5. Since for d > 1.5, we do not have an 
accurate approximation for the spectral function for the unparticle propagator, we technically 
have not ruled out the possibility that the effects of loop corrections to the propagator on 
g — 2 are large in this case. 
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Table 2: The relative difference 6(a^) between the unparticle contribution to 
muon g — 2 with and without vacuum polarization corrections to the unparticle 
propagator. An entry in the first column gives the choice of parameters used 
for that row, with ji and A3 given in GeV. The second column gives a^f from 
unparticles with no vacuum polarization corrections. The remaining column 
gives the fractional difference S(aM) — (aM — ott where a^f takes into 
account vacuum polarization correction to the unparticle propagator from all 
charged SM fermions. 
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Figure 5: The unparticle contribution to the muon anomalous magnetic moment 
(a") (left panel), and to the electron anomalous magnetic moment (a^) (right 
panel), as a function of the conformal breaking scale /z. Here, A3 = 1 TeV and 
we show curves for d — 1.1, d — 1.5, and d = 1.9. Shown for comparison in each 
case is the experimental limit (see text). 




(a) (b) 

Figure 6: Lower bounds on d and fi from muon (left panel) or electron (right 
panel) g — 2 for various values of A3. The region above the curves is allowed. 
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calculate the effects of this term we utilized an approximation based on disper- 
sion integral techniques, which is exact when there are no vacuum polarization 
corrections to the propagator. This approximation is accurate for d < 1.5; it 
breaks down for greater values of d. 

We found that adding the width-like term to the propagator significantly 
changed the the value of the propagator, but made little difference in the cal- 
culation of g — 2. Specifically, for the points we analyzed numerically using 
our dispersion integral approximation, we found that adding this term never 
changed the value of g — 2 by more than 1%. It would be interesting to learn 
whether this is the case for other observables or whether g — 2 is exceptionally 
insensitive to this term in the propagator. 

We also obtained limits on the unparticle dimension d and the conformal 
scale /i for various values of the unparticle scale A3 by assuming that the di- 
mensionless coupling c is equal to unity. In the muon case, we find that even 
for A3 as high as 1000 TeV, there are strong bounds on fi for d < 1.3. In 
the electron case, the bounds are much weaker, though we still find that for 
A 3 = 1000 TeV, fj, > 1 GeV for d < 1.2. 
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